The Diffuse Supernova Neutrino Background is detectable in Super-Kamiokande 
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The Diffuse Supernova Neutrino Background (DSNB) provides an immediate opportunity to study 
the emission of MeV thermal neutrinos from core-collapse supernovae. The DSNB is a powerful 
probe of stellar and neutrino physics, provided that the core-collapse rate is large enough and that 
its uncertainty is small enough. To assess the important physics enabled by the DSNB, we start 
with the cosmic star formation history of Hopkins & Beacom (2006) and confirm its normalization 
and evolution by cross-checks with the supernova rate, extragalactic background light, and stellar 
mass density. We find a sufficient core-collapse rate with small uncertainties that translate into a 
variation of ±40% in the DSNB event spectrum. Considering thermal neutrino spectra with effective 
temperatures between 4-6 MeV, the predicted DSNB is within a factor 4-2 below the upper limit 
obtained by Super-Kamiokande in 2003. Furthermore, detection prospects would be dramatically 
improved with a gadolinium-enhanced Super-Kamiokande: the backgrounds would be significantly 
reduced, the fluxes and uncertainties converge at the lower threshold energy, and the predicted event 
rate is 1.2-5.6 events yr _1 in the energy range 10-26 MeV. These results demonstrate the imminent 
detection of the DSNB by Super-Kamiokande and its exciting prospects for studying stellar and 
neutrino physics. 

PACS numbers: 97.60.Bw, 95.85.Ry, 98.70.Vc 



I. INTRODUCTION 



Core-collapse supernovae (CCSNe) occur at a rate of 
several per second in the Universe, each releasing a pro- 
lific ~ 10 58 neutrinos and antineutrinos. Their detec- 
tion provides a rich bounty for stellar and neutrino stud- 
ies, shown by the detection of the neutrinos from SN 
1987A 0, 0, I 1 1, S, 0, H, S M HI O- While a 
core-collapse supernova in the Milky Way would easily 
be detected in neutrinos, the occurrence rate is only < 
3 per century [ID, 0, EH • Proposed neutrino detectors 
should be able to detect supernovae up to 10 Mpc away 
with an occurrence rate of ~ 1 per year [l6l \vj\ . The 
vast majority of supernovae are therefore undetectable. 
However, the cumulative emission from all past core- 
collapse supernovae, which forms the Diffuse Supernova 
Neutrino Background (DSNB), h as p romising detection 
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130), |3i|, 122,123 1- The Super-Kamiokande (SK) limit of 
(t>{E Ve > 19.3 MeV) < 1.2 cm" 2 s" 1 || on the DSNB 
flux is already close to theoretical predictions [13, H2] • 

Predicting the DSNB for a given supernova neutrino 
emission model requires knowledge of the rate of core- 
collapse supernovae. In the past, this was not well known, 
and various studies provided insights on the supernova 
rate from the DSNB [H H3, HI- Indeed, the SK limit 
on the DSNB flux is strong enough to rule out some su- 
pernova rate evolution models, assuming a fiducial neu- 
trino emission model. However, our understanding of the 
cosmic star formation history (CSFH) has been greatly 
augmented by improved direct measurements in different 



wavebands and redshifts (see, e.g., Refs. [H, |37| and ref- 
erences therein). Cross-checks with other well-measured 
observables are now constraining, so that the CSFH is 
well determined by methods other than the DSNB. Thus 
it is both timely and important to study the prospects of 
the DSNB for probing stellar and neutrino physics. 

In this paper we start with the CSFH of Hopkins & 
Beacom (hereafter HB06 [13]) which is based on a com- 
pilation of recent data, and assess it by cross-checks, 
with the aim of evaluating the uncertainties that are car- 
ried forward into DSNB predictions through the core- 
collapse supernova rate. We henceforth refer to these as 
the astrophysical inputs and uncertainties on the DSNB. 
On the other hand, we refer to the supernova neutrino 
emission and neutrino properties as the emission inputs 
and uncertainties on the DSNB. 

As cross-check material we consider measurements of 
the rate of core-collapse supernovae, which have been sig- 
nificantly updated [38|, l39j, |4fJ, [4l|, |42| , the extragalactic 
background light, which records the total stellar emission 
over all time (for a recent review, see, e.g., Re f. 14311 ). and, 
finally, the stellar mass density (see e.g., Ref. [44j and ref- 
erences therein). While our approach is similar to previ- 
ous studies such as Ref. |30j , we perform novel checks and 
with higher precision. In particular, our analysis delivers 
fiducial inputs with unprecedentedly small uncertainties. 

Using our constrained astrophysical inputs, we find 
that the DSNB uncertainty is dominated by the emis- 
sion inputs, demonstrating the potential to study stel- 
lar and neutrino physics using the DSNB. Furthermore, 
taking into account both astrophysical and emission in- 
puts, we find that the predicted DSNB is within at 



most a factor ~ 4 of the SK limit set in 2003. For a 
6 MeV thermal spectrum, typical of scenarios with neu- 
trino mixing, this factor reduces to ~ 2. At the lower de- 
tection threshold energy for a gadolinium-enhanced SK 
(~ 10 MeV (45l H6[), the combined uncertainty on the 
predicted DSNB is remarkably only a factor 2, with a 
event rate of 1.2-5.6 events yr -1 in the energy range 
10-26 MeV. Thus, the core-collapse rate and the neu- 
trino emission per supernova are large enough to allow 
the imminent detection of the DSNB in SK. The suc- 
cessful detection will confirm the ubiquitous emission of 
neutrinos from core-collapse supernovae and initiate the 
much-anticipated study of stellar and neutrino physics, 
while a non-detection would require new stellar or neu- 
trino physics. 

The paper is organized as follows. In Sections [TT| and 
IIIII we discuss the astrophysical and emission inputs for 
the DSNB. We discuss the DSNB detectability and make 
future predictions for a gadolinium-enhanced SK in Sec- 
tion IIV1 and finish with conclusions in Section [V] We 
adopt the standard ACDM cosmology with Q, m = 0.3, 
Ov = 0.7, and H = 70 km s" 1 Mpc" 1 . 



II. ASTROPHYSICAL INPUTS 

For a given supernova neutrino spectrum, the key in- 
put for calculating the DSNB flux is the history of the 
rate of CCSNe (this includes Type II and the subdomi- 
nant Type Ib/c supernovae). The CCSN rate is directly 
related to the birth and death rate of massive stars. In 
recent years, data on the star formation rate (SFR) have 
improved both in breadth and sophistication, leading to 
an unprecedented understanding of the CSFH. In this 
section, we start with the latest CSFH and cross-check it 
with the measured CCSN rate, extragalactic background 
light (EBL), and stellar mass. Importantly, these checks 
probe different stellar life phases and are sensitive to dif- 
ferent stellar masses, making them complementary. Thus 
we are able to self-consistently assess the astrophysical 
inputs for the DSNB. 



A. Cosmic star formation history 

The SFR is most often derived from measurements of 
living massive stars. The measured luminosities, together 
with knowledge of their masses and lifetimes, gives their 
birth rates. Since the most massive stars have the short- 
est lifetimes, they provide a measure of the most recent 
star formation activity. In practice, the observed lumi- 
nosities are corrected for dust and the total SFR (over the 
entire stellar mass range) is derived by extrapolation to 
lower masses using the initial mass function (IMF). The 
calibration is done by use of a stellar population code 
that calculates the radiative output from a population of 
stars given an IMF. 




3 



redshift 



FIG. 1: Redshift evolution of the comoving SFR density. Data 
derived from various indicators are included as labeled, all 
scaled to a Salpeter IMF [54|. The majority are from the com- 
pilation of HB06 13711 . with additional data from LBG [551 ] and 
Ha measurements [5g| . We plot the fiducial CSFH (thick solid 
curve) and our generous adopted uncertainty range (thin solid 
curves). The curves take into account SFR data at higher red- 
shifts 37, 57] that are not shown here, where we focus on the 
lower redshifts most relevant for the DSNB. 



We compute calibration factors using the PEGASE.2 
stellar population code [47j . which contains a careful 
treatment of stellar physics. We assume constant SFR 
bursts of 10 s years, a close binary fraction of 0.05, evo- 
lutionary tracks with stellar winds, the supernova model 
B of Woosley & Weaver [48} , and a constant metallicity 
of Z = 0.02 (i.e., solar). The parameter that the results 
are the most sensitive to is the star formation duration. 
For integrated measurements of galaxies, it is usually 
appropriate to assume the SFR has remained constant 
over time scales that are long compared to the lifetimes 
of the dominant UV emitting population (~ 10 8 years) 
[49l ]; however, there are calibration uncertainties of a few 
tens of percent (see Table [I]). In contrast, studies suggest 
that the mean metallicity of star formation gas is close 
to the solar value for redshifts of a few [5(| HlJ , resulting 
in calibration uncertainties of ten percent or less. The 
other parameters yield variations of only a few percent. 
We note the agreement of PEGASE.2 outputs with other 
stellar population synthesis codes demonstrated in e.g., 
Refs. (5|, Hi]. 

We calculate results for three IMFs: the traditional 
steeper Salpeter IMF [54[, an intermediate Kroupa IMF 
(HI, and a shallower Baldry-Glazebrook (BG) IMF 59]. 
We define the IMF as ip(M) = dN/AM so that ip(M)dM 
gives the number of stars in the mass range M to M + 
dM. The slopes ip(M) oc are shown in Table [0 

We take the SFR compilation and fit of HB06 as our 
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TABLE I: Calibration factors, /bed = Lhoi/p* and /uv = 
Lvv[p*, calculated using PEGASE.2 [H, for three IMFs [HI, 
[H^,[53] and three epochs (yr). For a particular IMF, the star 
formation duration causes an uncertainty of a few tens of 
percent about the central 10 s yr value. This contributes to the 
scatter in the SFR for a given IMF. Note /uv is determined 
at A = 0.2 (jm since the UV output is nearly constant over 
wavelength for the durations considered. 



TABLE II: CSFH parametric fit to the form of Eq. © 
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starting point and add recent data. The HB06 data con- 
sist of various SFR indicators, including UV measure- 
ments from SDSS [H, GALEX [HI [pi and COMB017 
[63| . far-infrared (FIR) from Spitzer 64], and high red- 
shift measurements from the Hubble Ultra Deep Field 
(UDF) [65[; to this we add recently-derived data from 
Lyman Break Galaxies (LBG) [H| and from Ha emission 
[5a] , These are shown in Fig. [T] as a function of redshift, 
all scaled to a Salpeter IMF. Assuming the Kroupa or 
BG IMF results in values that are lower by an overall 
factor ~ 0.66 and ~ 0.55, respectively. 

Fig. Q] demonstrates the overall consistency between 
SFR densities estimated by different indicators, over a 
wide range of redshifts. In general, the systematic scat- 
ter due to different indicators dominates over the for- 
mal uncertainties from the calibrations of each indicator. 
However, at low redshifts (z < 1), the scatter decreases 
and approaches the calibration uncertainties. To further 
constrain the CSFH in the future would require calibra- 
tion uncertainties to be examined. We note that red- 
shift dependent dust corrections have been applied. At 
z < 1, the UV and FIR measurements are combined, 
while for 1 < z < 3, a constant dust correction is made; 
see Ref. [36[ for further details. 

To compare the CSFH to other observables, it is useful 
to define an analytic fit. We adopt a continuous broken 
power-law as in Ref. [57j , 



P*{z) = Po 



(1 + *) 



l + z 
B 



1 



C 



1/7) 



(1) 

where po is the normalization, B and C encode the red- 
shift breaks, the transitions are smoothed by the choice 
rj ~ —10, and a, p and 7 are the logarithmic slopes of 
the low, intermediate, and high redshift regimes, respec- 
tively. The constants B and C are denned as 

B = (1 + zx) 1 -^, (2) 
C = (l + z 1 )(^ tt )^(l + z 2 ) 1 -^, (3) 
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"Shown for the Salpeter IMF. For the Kroupa and BG IMFs the 
normalization po decreases by a factor ~ 0.66 and ~ 0.55, respec- 
tively; the overall shape is not greatly affected (see, e.g., Table 2 of 
Ref. [|3l). Units of p are in M yr _1 Mpc -3 . 



where z\ and Z2 are the redshift breaks. We adopt the 
fiducial CSFH fit obtained by combining the HB06 com- 
pilation with new measurements at high redshift derived 
from gamma-ray bursts [H3, [fiH] (note that high-redshift 
data are not shown in Fig. [I} . We define a generous en- 
velope that takes into account the scatter in the data, as 
shown in Fig. [T] The parameters are given in Tabic [TTJ 



B. Rate of core-collapse supernovae 

The number of stars per unit mass undergoing core 
collapse is dependent on the IMF and mass range of stars 
that lead to core collapse. The predicted comoving CCSN 
rate history, i?ccSN(z), follows directly from the CSFH, 



-RccsnO) = p*(z) 



J*°HM)dM 
j^ MiP(M)dM' 



(4) 



where the ratio of integrals is 0.0070/M Q , O.O1O9/M , 
and 0.0132/M Q for the Salpeter, Kroupa, and BG IMFs, 
respectively. 

The lower mass threshold is the most important pa- 
rameter for calculating i?ccsN- In general, it is difficult 
to predict accurately from theory, because stellar prop- 
erties change rapidly between ~ 6-10 M Q ; it is usually 
assumed to be 8 M Q . On the other hand, the mass can 
be determined from direct identifications of progenitor 
stars from pre-explosion imaging. The most recent study 
places the mass threshold for Type II-P supernovae at 
8.5^1 5 M Q 41]. The uncertainties would translate to an 
uncertainty of about ten percent in i?ccSN- However, one 
should keep in mind that explosions of O-Ne-Mg cores 
may result from a lower mass range that leads to a dif- 
ferent variety of CCSN. The upper mass limit is less im- 
portant because of the steep IMF, and only affects the 
results at the percent level. Theoretically, prompt black 
hole formation occurs above some critical mass; the pre- 
cise value depends strongly on properties such as rotation 
and metallicity, but is thought to lie around 30-50 M Q 
[67l ] . Such explosions will generally not produce optical 
supernovae but will likely produce neutrinos. Observa- 
tionally, a progenitor mass of at least 40 M has been 
claimed to have been inferred from a neutron star in a 
very young cluster f68j . 
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FIG. 2: Evolution of the core-collapse supernova rate (top 
panel) and Type la to CCSN ratio (bottom panel), shown 
against distance in Mpc (for z < 0.1) and redshift (for 
z > 0.1). Data as labeled; filled symbols indicate data or 
limits we adopt, and empty points indicate published rates 
we treat as lower limits. Error bars show statistical errors, 
except for horizontal bars on Dahlen et al. points that show 
bin size. In the top panel, predictions from the CSFH are 
shown by the solid curve, where the uncertainty band propa- 
gates from Fig. [Tj The SFR measurements, scaled to a CCSN 
rate, are shown in the background (small light brown squares) . 
The BG IMF has been adopted; using other IMFs causes 
differences of only a few percent. In the bottom panel, the 
dashed lines denote conservative upper limits from the non- 
observation of any Type la supernovae within 10 Mpc [l7| . 
The 10 Mpc point has been very conservatively selected from 
Ref. [l7( ] . and the rate has been further corrected downwards 
(see text). 



It should be noted that the predicted i?ccSN is largely 
insensitive to the IMF. This is because the effects of the 
IMF nearly cancel out between p* and the integrals in 
Eq. iQJ. For instance, a shallower IMF contains rela- 
tively more massive stars but has a smaller p*. Numer- 
ically, e.g., i?ccSN ~ 0-96-Rccsn- This is a natural out- 
come because both the CCSN and the SFR indicators 
are associated with the same massive stars. Thus, the 
CCSN rate allows an almost IMF-independent check of 
the CSFH. 

In the top panel of Fig. [5] we show the predicted CCSN 
rate, using the BG IMF, for the fiducial CSFH (thick 
solid curve) as well as the upper and lower CSFH (thin 
solid curves). Using the Kroupa or Salpeter IMF pro- 
duces variations of only a few percent, as expected. For 
data, we plot rates derived from supernova compilations 
[ftL Pill ], extended supernova surveys using the Hubble 
Space Telescope [42J, searches in nearby galaxies 
[3 and the Southern inTermediate Redshift ESO 
Supernova Search (STRESS) 0. We also show in the 



background the SFR measurements scaled via Eq. ((4]) to 
a CCSN rate using the BG IMF. 

For the local Universe, CCSN rates are derived from 
compilations of supernovae occurring within 10 Mpc over 
a time period of 10 years I17| . and within 28 Mpc over 
a time period of 10.5 years [41(. These compilations in- 
clude Types la, II, lb, and Ic supernovae, the latter three 
of which are used to derive the CCSN rate. Due to the 
known incompleteness of the compilations, both should 
be treated as lower limits on the CCSN rate. In ad- 
dition, for the 10 Mpc point, we conservatively exclude 
sources that might inadvertendly be outside 10 Mpc and 
we further exclude transients of peculiar origin (e.g., SN 
2002kg, SN 2008S and the 2008 transient in NGC 300) 
[70T ]. which reduces the number of CCSN to 13. Also, the 
10 Mpc point has been corrected downwards for the lo- 
cal increase in SFR density. We apply a simple correction 
assuming the SFR density traces the galaxy density; the 
galaxy over-density in 10 Mp c is about a factor 2 accord- 
ing to galaxy catalogues [7l| (see also Fig. 2 of Ref. [72j]). 
The correction does not affect the Type la to CCSN ratio 
that we discuss below. 

Since the lifetimes of massive stars are short, the shape 
of the CCSN rate evolution must follow that of the 
CSFH. On the other hand, the normalization needs to 
be checked. We find excellent agreement between predic- 
tions from the CSFH and the most recent data from the 
survey of Dahlen et al. using the Hubble Space Telescope 
[42l |. These data have been derived by periodically ob- 
serving the same patch of sky, locating supernovae within 
a volume only limited by flux. The new data (red solid 
squares) [42| are updates of their previous data (blue 
empty circles) (39j , and have much better statistics and a 
more detailed analysis (we plot the data which have been 
corrected for dust). We find further agreement with re- 
cently published supernova compilations within 10 Mpc 
[13] and 28 Mpc [4l| , even though they are certainly in- 
complete, especially in the Southern sky, confirming the 
minimal normalization over a wide distance scale. 

On the other hand, data from Cappellaro et al. [H, [69| 
and STRESS [H fall short of the trend of other data. 
Unlike in Dahlen et al., these surveys search for super- 
novae by periodically observing a pre-selected sample of 
galaxies in a given field. While a large number of galax- 
ies are selected, even this is likely incomplete, as small 
galaxies are often under-sampled. In addition, as the 
authors clarify, host galaxy extinction is an important 
uncertainty on supernova surveys [4rJ |. The surveys of 
Cappellaro et al. do not include host galaxy extinction, 
and as the authors state, they most likely have too few 
CCSN [H. 

The ratio of Type la to CCSN also assists in assessing 
the results of supernova surveys. Due to their brighter 
nature, Type la supernovae are more easily detected, so 
that the ratio of Type la to CCSN will increase if CCSN 
are increasingly missed. In the bottom panel of Fig. 
we show the ratio of raw supernova counts for Ref. [4l|, 
and the ratio of reported Type la and CCSN rates for 



5 



other studies. From the non-observation of any Type la 
supernovae within 10 Mpc - of the supernovae consid- 
ered in the 10 Mpc compilation, none were Type la [13 
- we determine 90% and 99% upper limits on the ratio, 
shown by dashed lines. We apply these upper limits over 
all distances and redshifts; although the ratio in princi- 
ple evolves with time due to Type la time-delay effects, 
the evolution is not strong for global samples and z < 1 
(see, e.g., Refs. [zl, [z3)' an d moreover, a delayed com- 
ponent to the la rate relative to the CCSN rate would 
only make the limits stronger by increasing the local ra- 
tio. We therefore use the dashed lines as a conservative 
indicator for the true ratio, for simplicity. We see from 
Fig. [2] that more CCSN are being missed in the data 
we treat as lower limits (empty symbols). For example, 
the large ratio of Ref. [69( predicts ~ 5 Type la within 
10 Mpc, in disagreement with observations [17] . Alter- 
natively, taking the 90% upper limit as described above 
and respecting the reported Type la supernova rates, the 
CCSN rates (empty symbols) increase to values more in 
line with theoretical predictions from the CSFH. 

We therefore conclude the data confirm the CSFH nor- 
malization over a large range of distances. Moreover, the 
range is precisely that of interest for the DSNB (almost 
all of the detectable flux comes from z < 1). 

As an aside, while predictions include all core collapses, 
including those with little or no optical signals due to 
the prompt formation of a black hole or due to dust ob- 
scuration, the data are derived from optical supernovae 
only. The excellent agreement between prediction and 
data suggests that dark core collapses comprise a mi- 
nority of all core collapses. Alternatively, the fraction 
of optically-dark but neutrino-bright collapses can be in- 
creased by adjusting the CSFH and mass range for core 
collapse; however, this would be accompanied by a cor- 
respondingly larger DSNB flux, which is constrained by 
the tight experimental limits discussed below. 

We note in this context the red supergiant problem 
noted by Smartt et al. I4l| . They find a shortage of 
Type ITP supernovae associated with red supergiants in 
the mass range 17-25 M Q , despite clear evidence that 
these massive supergiants exist. These stars constitute 
~ 10% of all massive stars leading to core collapse, and 
although slightly smaller in mass than theoretical pre- 
dictions, they could be dark core collapses. Monitoring 
nearby red supergiants can reveal the occurrence rate of 
dark core collapses, as proposed by 



C. Extragalactic background light 

It is widely accepted that the EBL is the record of the 
total stellar emission over all time. The observed EBL is 
dominated by two peaks of comparable energy density. 
The first peak in the optical to near infrared (NIR) is 
powered by direct starlight, while the second in the FIR 
is dominated by starlight that is absorbed and re-emitted 
by dust. Since circumstellar dust absorbs most efficiently 



in the UV, the power source of the FIR peak is primarily 
massive stars. The relative energy densities of the two 
peaks is a testament of the heavy dustiness of the present 
and past Universe (for a review, see e.g., Ref. |43j). 

We check the CSFH by using the total observed EBL. 
In principle, if one had precise knowledge of dust - 
its quantity, properties and evolution - one could use 
the spectral shape of the EBL to probe different stellar 
mass regimes. However, a full treatment is beyond the 
scope of this paper and we simply exploit the total EBL 
as a calorimeter of the total energy radiated by stars. 
Although this includes contributions from less massive 
stars, half of the total EBL is powered by stars with 
masses M > 3M for the BG IMF. 

Direct measurements of the EBL in the optical to NIR 
are complicated by foreground contamination [43j . In- 
direct limits have been placed from integrated galaxy 
counts and from the opacity of the Universe to gamma 
rays using distant TeV blazars. We caution however that 
gamma-ray constraints depend on the assumed intrinsic 
blazar spectral index, so the result may be weakened in 
light of unknown acceleration mechanisms [76l WK . [78l . 
[79l . l80l | . In Fig. [3] we present the latest data on the 
EBL, including direct measurements (empty symbols) 
[HI m, HI Hi IM HE [13j8ij8ij[ower limits from galaxy 
counts (filled symbols) |fl,|9l|,S M, M, S M, [Hf , 
and up per limi t s fr om gamma-ray attenuation (thick 
curves) pi flOOL fToij . 

The systematic uncertainties on galaxy photometry 
and zodiacal light subtraction suggests that a curve be- 
tween the measurements and counts represents the most 
appropriate assumption for the total observed EBL. We 
show three EBL regions: the minimum, which essen- 
tially traces the galaxy counts, the nominal, which re- 
spects the gamma-ray constraints, and the maximum, 
which lies within most of the data error bars. We do not 
consider the 1-4 micron diffuse sky emission detected by 
IRTS [l03| to be of extragalactic origin, since it is most 
likely caus ed by emission from interplanetary dust par- 
ticles 105]. The total EBL from these regions are 52, 
73 and 99, all in units of nWm~ 2 sr . Ou r est imates 
are slightly higher than those of Refs. [9ll 106], who 
find a best estimate of 60 nWm~ 2 sr _1 . On the other 
hand, they ar e mor e in line with recent estimates of 60-93 
nWm^sr" 1 flOp, 45-170 nWm" 2 sr" 1 and 50-129 
nWm -2 sr -1 [531 ]. 

To calculate the EBL for a given CSFH, we need the 
spectral luminosity density e{y, z) as a function of z, mea- 
sured since the epoch z* when stars first turned on. As- 
suming the EBL arises dominantly from stellar radiation, 
e is given by jl08l | 



e(y,z) 



p*(t)dt / L(u,M,t/)ip(M)dM, (5) 



where L(v, M, t') is the luminosity per unit mass of a star 
of initial main-sequence mass M at time t' = t z —t, ip(M) 
is the IMF in the range 0.1-100 M and M(t') is the 
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FIG. 3: The observed EBL spectrum. Filled symbols are 
based on integrated galaxy counts, while empty symbols and 
the black lines at high wavelength represent absolute measure- 
ments. In the UV to NIR, we also show the gamma-ray upper 
limits placed by H.E.S.S. (red solid curve) MAGIC (blue 
dashed curve) [lOll l , and Edelstein et al. 2000 (filled downward 
triangle at 0.1 /jm) |102| ]. We show three shaded regions repre- 
senting the minimum (dark shading) , nominal (light shading) 
and maximum (lightest shading) EBL. The integrated total 
EBL from these are listed at the top-right. 
The integrated galaxy counts shown include (in order of in- 
creasing A) Gardner et al. 2000 (red cicle) [9(j, Madau & 
Pozzetti 2000 (green square) Levenson & Wright 2008 

(blue diamond) [9^], Fazio et al. 2004 (yellow right-pointed 
triangle) [93|. Metcalfe et al. 2003 (red up-pointed triangle) 
[||, Elbaz et al. 2002 (blue down-pointed triangle) [H , Pa- 
povich et al. 2004 (green right-pointed triangle) [9q |. Dole et 
al. 2006 (yellow circle) and Frayer et al. 2006 (purple 
square) 98] . Absolute measurements shown include (in order 
of increasing A) Cambresy et al. 2001 (dark blue right-pointed 
triangle) j85|], Wright 2001 (purple left-pointed triangle) [§f|, 
Dwek & Arendt 1998 (dark green square) [87j], Lagache et 
al. 2000 (green diamond) Hf, and Hauser et al. 1998 (red 
cirlce) [gg|. For the results of Bernstein et al. [8ll. l82l, l83l [8^ . 
we show 2a upper limits (filled downward-pointing trian- 
gles), nominal values (empty circles), and lower limits (filled 
upward-pointing triangles) , all connected by dashed lines. We 
show the diffuse sky m easu rements of Matsumoto et al. 2005 
(small empty circles) [103| ]. and the FIR measurements by 
Fixsen et al. 1998 (filled black circles) [104| ] connected by black 
lines. 



initial main-sequence mass of a star with a lifetime of t'. 
We calculate this quantity using the PEGASE.2 stellar 
population code, with the same assumptions as described 
in Section [II Al but with evolution of metallicity, with an 
initial value of Z = 0.001 at z = 10. The specific EBL 
intensity at the observed frequency v is then the 
integral of e{y' , z), from the comoving volume element at 



z, over redshifts, 



1(1/) 



c 

4tt 



e(i/,z) 



At 



dz, 



(G) 



where v' = v{l + z) is the frequency at emission and 
|dz/di| = H (l + z)[n m (l + zf + 0^/2. 
The values of the calculated total EBL are 



j+31 



951^ nWm^sr" 1 , 88lf 8 nWm^sr" 1 , and 78l 24 
nWnr 2 sr _1 for the Salpeter, Kroupa, and BG IMFs, 
respectively. While the CSFH differs by almost a factor 
2 between the Salpeter and BG IMFs, the calculated to- 
tal EBL differs by much less, because the total EBL is 
dominated by relatively higher mass stars. On the other 
hand, steeper IMFs have more low mass stars which live 
long and pile up, which works to increase the total EBL. 
These dependencies on the IMF have in fact been studied 
by various authors to constrain the IMF (e.g., Ref. ■53]). 
Of importance to us is the agreement between observa- 
tions and predictions, when calculated using recent IMFs 
with shallower slopes and suppression at the lower mass 
end (BG IMF). 

We note that this result contrasts with a recent study 
which found predictions that were smaller than observa- 
tions [53]. In their study, the predicted EBL was ~ 50 
nWm _2 sr _1 for the BG IMF and they estimated the 
observed EBL to be 77 nWm _2 sr _1 . Our study differs 
in three aspects. First, by including the gamma-ray con- 
straints, the updated observed EBL is slightly smaller. 
Second, the updated CSFH is somewhat larger. Finally, 
we include consistent evolution of metallicity in our calcu- 
lations, whereas the authors in (53| assumed a constant 
Z = 0.02 over all redshifts. Lower metallicity leads to 
less mass loss and hence higher time-integrated radiative 
output. These three factors result in us obtaining better 
agreement. 



D. Total stellar mass 

Integration of the CSFH over redshift with appropri- 
ate corrections for stellar mass loss yields the stellar mass 
density. This quantity can be independently measured 
using galaxy surveys, which are often coupled to NIR ob- 
servations as a proxy for stellar mass. Therefore, it also 
provides an independent check of the CSFH. The com- 
parison has the property of probing a lower stellar mass 
range than the EBL. Numerous studies ha ve been made , 
with results varying from good agreemen t [l09l IllOl llll| 
to the CSFH over-producing stars [3?], Ill2ll. In these 
comparisons, the IMF plays a critical role jl!3| . A recent 
detailed study shows that the CSFH and observations of 
stellar mass density are in good agreement for redshifts 
z < 0.7 (44j. In a subsequent paper, the authors find that 
if the IMF is constant in time, the best-fit IMF slope is 
2.15 [ll4| . which is the same as our adopted BG IMF. 
Although the check becomes complicated by large scat- 
ter in measurements at redshifts above 0.7, the studies 
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illustrate the overall consistency of the CSFH and stellar 
mass in the redshift of our interest and the preference for 
a shallow IMF. 



E. Consistency of the CSFH and observations 

Besides systematic uncertainties that contribute to the 
scatter in the CSFH, another issue is the CSFH normal- 
ization uncertainty arising from dust correction. As we 
discuss below, the true CSFH cannot be smaller by a 
factor at most ~ 2, and even this seems unlikely. 

Various studies have shown that the EBL is dom- 
inated by stars, with little contribution from non- 
nucleosynthesis energy sources such as active galactic nu- 
clei I53l. Ill5l. Hence, we infer the minimum CSFH from 



the minimum observed EBL, « 50nWm sr . Requir- 
ing the minimum observed EBL to be explained, the true 
CSFH could be at most a factor 95/50 ~ 2 smaller than 
the fiducial CSFH. 

However, this would require several unlikely changes. 
First, the Salpeter IMF is disfavored by other observ- 
ables such as the stellar mass density. Second, the dust 
correction applied to the CSFH is typically a factor 2- 
3. Invoking a true dust-correction of 1-1.5 would imply 
almost negligible true dust, in conflict with the observed 
FIR EBL peak. Third, to obtain consistency with the 
CCSN rates, one would need to either increase the mass 
range for CCSN or invoke almost negligible true dust 
correction, and neither seems plausible. 

Therefore, we conclude that the CSFH cannot be 
smaller by even a factor ~ 2; the maximal reduction fac- 
tor is perhaps 78/50 ~ 1.5. This sets the lower limit of 
the CSFH normalization. We remind the reader that the 
DSNB detectability is not directly affected by the value 
of this factor, since the flux is directly normalized by the 
CCSN rate and is independent of uncertainties associated 
with low-mass star formation. 



shifted, over cosmic time [ll6j . 



dz, (7) 



where E' = E(l + z). 

Progenitors over a wide range of masses lead to sim- 
ilar neutron star masses and hence neutrino emissions 



1171 ] . The dominant neutrino emission occurs during the 



Kelvin-Hclmholtz cooling phase, when the newly formed 
hot a nd dense protoneutron s tar c ools to a neutron star 
[l 18l | (see also Refs. [l!9| | and [l20l |). Neutrinos and anti- 
neutrinos of all flavors are produced (y e , D ei and v x \ 
where v x refers to v^, v T , and their antiparticles) , and 
each species carries away an approximately equal fraction 
of the total energy, i?* ot ss 3 x 10 53 erg. Their spectra are 
to a good approximation thermal; we summarize some 
temperatures from numerical supernova simulations in 
Table IIIII (where we use i to explicitly denoted quan- 
tities at production). The hierarchy T* e < < 
reflects the different radii at which each neutrino species 
decouples, which in turn arises from the relevant neutrino 
interactions. 

The observed v e spectra outside the star are linear 
combinations of the neutrino spectra at production, ow- 
ing to neutrino mixing. The time-integrated v e spectrum 
per supernova is well approximated by t he Fermi -Dirac 
distribution with zero chemical potential [l 19L Il20j , 



dN, 



dEi-, 



El ot 120 E, 



Us. ( e E k 



/Tp, 



(8) 



where T Pc is the effective D e temperature outside the star 
after neutrino mixing. This is the temperature that is 
measured by neutrino detectors and we therefore use it 
for predictions. The effective temperature contains in- 
formation on stellar and neutrino physics, and it is a 
separate problem to work backwards from the effective 
spectrum to the initial spectra, taking into account the 
effects of neutrino mixing. 



III. DSNB PREDICTIONS 

In the previous section we explored the consistent pic- 
ture of stellar birth, life, and death. As a result, we ob- 
tained the fiducial astrophysical input and uncertainties. 
In this section, we introduce the DSNB and discuss the 
neutrino emission per supernova and neutrino properties, 
i.e., the emission inputs and uncertainties. We then make 
predictions for the DSNB flux. 



A. DSNB formalism 

The predicted DSNB number flux, over Air, is calcu- 
lated by integrating i?ccSN(^) multiplied by the neu- 
trino emission per supernova, dN/dE, appropriately red- 



B. Neutrino emission per supernova 

It is important to address the range of expected 
neutrino emission per supernova, because we integrate 
over the entire CCSN population. Potential processes 
that affe c t th e neu t rino emission are both microphys- 
ical [ml, EH HH EH and macrophysical EH EM 
Il27l Il28l | . Ne utrino mixing also plays an important role 
[28l . I129L Il30f ] . Although a body of predictions is given 
in the literature, it should be emphasized that numer- 
ical simulations usually do not explode, and that even 
the most state-of-the-art simulations do not reach be- 
yond the accretion phase at a few hundred milliseconds 
after bounce. We therefore try to be as general as possi- 
ble and focus on the time-integrated emission needed for 
the DSNB. 
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TABLE III: Flavor-dependent temperatures from some exam- 
ples of numerical supernova simulations in the literature. 



Author 


time" 


rpi b 


rpi b 


rpi b 


Ref. 


Myra & Burrows (1990) 


0.2 


3.3 


4.0 


8.0 


[121] 


Totani et al. (1998) 


0.5 


3.9 


4.9 


6.3 


[132] 




10 


3.5 


6.3 


7.9 




Rampp & Janka (2000) 


0.5 


2.3 


3.4 




[133] 


Liebendoerfer et al. (2001) 


0.5 


4.2 


4.6 


5.3 


[134] 


Mezzacappa et al. (2001) 


0.5 


3.6 


4.2 


5.2 


[135] 


Keil et al. (2003) 




3.7 


4.0 


5.2 


[136] 


Thompson et al. (2003) 


0.2 


2.9 


3.5 


4.5 


[137] 


Liebendoerfer et al. (2004) 


0.5 


3.7 


4.2 


4.5 


[138] 




0.5 C 


4.5 


4.8 


7.5 




Sumiyoshi et al. (2007) 


1.4° 


6.6 


7.0 


10 


[127] 



"postbounce time in seconds 
'temperatures at production in MeV 
c leads to black hole formation 



1. Total neutrino energies 

The total energy budget in all flavors of neutrinos is 
dictated by the binding energy of the final remnant, 

(9) 

where Mms and Rns are the neutron star mass and ra- 
dius. The neutron star mass is best estimated by mea- 
surements in high-mass X-ray binaries and binary pulsars 
containing a radio pulsar and a neutron star companion. 
The mos t likely value is « 1.4 M Q , with a range 1.2-1.6 
M [l39j | . The radii are more difficult to measure, but are 
consistent with being ~ 10 km. Note the higher masses 
measured in low-mass X-ray binaries are attributed to 
their longer mass accretion histories. 

The P e emission depends on how E v ot is partitioned be- 
tween neutrino flavors. In the numerical supernova sim- 
ulations of the Lawrence Livermore (LL) group, which 
successfully followed the simulation into the Kelvin- 
Helmholtz cooling phase (~ 18 s) [132], almost exact 
luminosity equipartition is seen throughout. The time- 
integrated total v e energy is 4.7 x 10 52 erg, in good agree- 
ment with E v ot /6. 

While the LL group succeeded in obtaining an explo- 
sion, it has been realized that they lacked neutrino pro- 
cesses now recognized as important. A review of the lit- 
erature shows that luminosity equipartition is not uni- 
versal: while L l v w Lp e appears to be robust, L v ^ varies 
by a factor 2-3 in either direction, depending on the evo- 
lutionary phase and on numerical methods. However, we 
remind the reader that these simulations only reach a few 
hundred milliseconds at most, when the bulk of neutrinos 
have not yet been emitted, and therefore say little about 
energy equipartition. 

In addition, the observed v e emission is a linear com- 



bination of v e and v x at production, due to neutrino 
mixing. Mixing scenarios, which are dependent on the 
neutrino mass hierarchy and oscillation parameters as 
well as the neutrino and stellar densities through which 
the neutrino s propagate, have been studied systemati- 
cally [H, HE [130, Eiil ■ These give the effective 
v e temperature a value in between the extremes corre- 
sponding to a pure v e at production and a pure v x at 
production. 

Observationally, analysis of D e 's from SN 1987A show 
that the observed D e energy budget is 3-6 x 10 52 erg 
[l42l Il43l ] , confirming the approximate energy equiparti- 
tion after neutrino mixing. Variations away from energy 
equipartition at production are likely reduced by the ef- 
fects of neutrino mixing as described above. 

In some cases, the total energy budget may be smaller 
than Eq. ([9]). One example is the effect of the equation of 
state for dense matter. Another possibility would seem 
to be the case of prompt black hole formation, where 
the neutr ino signal is expected to be abruptly cut short 
[l3cj Il44| . However, recent simulations show that black 
hole format ion is preceded by an increase in the neutrino 
luminosity [HI EH EH , as the gravitational energy of 
the rapidly accreting matter is released. The time inte- 
grated neutrino luminosity can be as high or higher than 
in normal neutron star formation. Another potential case 
is rotation, which could decrease Eq. © by bloating the 
neutron star, although the effects are most likely not sig- 
nificantly large for an integrated population (see the light 
curves converging with time in Fig. 2 of Ref. [1251 ]). These 
are all interesting physics to study using supernova neu- 
trinos. 



2. Average neutrino energies 

In general, the neutrino temperature is only weakly 
dependent on the neutron star parameters. As- 
suming a thermal spectrum, the neutrino luminos- 
ity is L v oc i? 2 X^. Since the total energy bud- 
get is well defined, it is useful to consider i?* ot w 
L u At, where At « 10s (R NS /10km) 2 (p/2p nm ) 2 / 3 is 
the time scale for the neutrino cooling phase [l46l . 
Il47| and p nm is the nuclear density. This gives 
T v ~ 5(M A rs/1.4M ) 1 / 3 (i? A r S /lOkm)- 3 / 4 MeV. Al- 
though this is a simplistic analysis, it shows that T„ is 
weakly dependent on the neutron star mass and radius. 
For example, the neutrino emission from a more massive 
neutron star-neutron star mer ger s ystem is comparable 
to that of a normal supernova [1481 ]. 

Indirect constraints on the v l x spectrum have been 
placed from the chemical abundance of ^-process in- 
duced elements. Since the inelastic v x scattering off nu- 
clei uniquely leads to the production of specific isotopes 
[l49j ] , the required v x flux can be inferred. For E v — 4- 
5.8 x 10 52 erg, T % v is constrained to lie between 4-7 MeV 

[Mfi5il. 

The DSNB v e spectrum can be calculated from the ob- 
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FIG. 4: DSNB flux spectrum for emitted neutrino spectra 
as labeled. For each spectrum, two curves are plotted rep- 
resenting the full range of uncertainties due to astrophysical 
inputs (the fiducial prediction lies in between). The shadings 
indicate backgrounds, with origins as labeled. Decays of in- 
visible muons and spallation products would be reduced in a 
gadolinium-enhanced SK, opening the energy region 10 MeV 
and above to a rate-limited DSNB search; see Fig. [5] 

served SN 1987A v e dlMIH. The SN 1987A spectrum 
constructed by fitting [3l| or by nonparametric methods 
[32| shows similar distortions from a thermal spectrum. 
For the purposes of DSNB detection, the spectrum near 
detector threshold is important, which as we discuss be- 
low is 20 MeV for the current SK and « 10 MeV in 
a gadolinium-enhanced SK. The spectrum of Ref. [32j 
is similar to a 5 MeV thermal spectrum for high ener- 
gies (E u > 30 MeV), and closer to a 4 MeV thermal 
spectrum at somewhat smaller energies (E v < 30 MeV). 
Therefore, the SN 1987A spectrum could be treated as 
similar to a 4 MeV thermal spectrum, which is on the low 
end of the theory predictions. The limited sample of the 
SN 1987A data means this spectrum serves as a guid- 
ance rather than the definitive neutrino spectrum per 
supernova; there is no good reason to invoke a smaller 
spectrum. 



C. DSNB flux prediction 

In Fig. |4] we show the predicted DSNB flux for a selec- 
tion of neutrino spectra as labeled. For each, the two 
curves correspond to the upper and lower i?ccSN in- 
puts from Fig. [2j The fiducial astrophysical input lies 
in between. The figure illustrates the relative sizes of 
the DSNB uncertainties originating from the astrophysics 
and emission inputs. The shading represents the relevant 




E e [MeV] 

FIG. 5: DSNB event rates at SK (flux spectra weighted with 
the detection cross section) against positron energy. Note the 
linear axis. We hatch in the 2003 upper limit by the Super- 
Kamiokande Collaboration, < 2 events (22.5 kton yr)~ in the 
energy range 18-26 MeV. The limit applies to all spectra (see 
text). In a gadolinium-enhanced SK, decays of invisible muon 
and spallation products would be reduced, opening up the 
energy range > 10 MeV for DSNB search (unshaded region). 



backgrounds which detections of the DSNB must com- 
pete against. We discuss in Sec. IIVCI how these back- 
grounds would be rejected in a gadolinium-enhanced SK. 

We adopt thermal spectra given by Eq. {5]) with effec- 
tive temperatures of 4, 6, and 8 MeV. These are effective 
temperatures, reflecting the range of initial temperatures 
shown in Table IIII1 as well as the range of neutrino mix- 
ing effects described in section IIIIBI We also show the 
neutrino spectrum directly reconstructed from the SN 
1987A neutrino data [Hj]. 

We see that with the current astrophysical uncertain- 
ties, the emission uncertainties dominate at E v > 20 
MeV. The DSNB in this energy range therefore provides 
information on the emission inputs. At lower energy the 
fluxes and uncertainties converge. 



IV. DSNB CONSTRAINTS AND DETECTION 
A. Super-Kamiokande constraints 

Super-Kamiokande is a water Cerenkov neutrino detec- 
tor with a fiducial volume of 22.5 kton, sensitive to DSNB 
through the inverse beta reaction D e p — ► e + n on free pro- 
tons (hydrogen nuclei). The positron energy faithfully 
represents the neutrino energy, and E e ~ E„ c — 1.3 MeV 
to zeroth order in M„. In 2003, the Super-Kamiokande 
Collaboration reported a stringent upper limit on the 
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integrated DSNB v e signal above E e — 18 MeV, using 
1496 days (4.1 yr) of data. The dominant background 
in this energy range was the decay of invisible muons, 
the fixed spectrum of which rises steeply with positron 
energy [26|, [34[ . The best signal to noise ratio is there- 
fore obtained at the lowest energy bin, 18-26 MeV, where 
the background rate is some ~ 3 events yr" 1 . By search- 
ing for DSNB neutrinos over detector backgrounds, the 
Super-Kamiokande Collaboration has limited the DSNB 
flux to be (f>(E Pe > 19.3 MeV) < 1.2 cm" 2 s" 1 . In a more 
physically relevant form, the limit is < 2 events (22.5 
kton yr)" 1 in the lowest energy bins [34. |46|. 



B. DSNB event spectrum 

In Fig. [5] we show the predicted event rates (flux spec - 
trum weighted with the detection cross section |l52l . ll53| ) 
as a function of positron energy, which is the observed 
quantity in SK. We account for corrections of order 1/M p 
in both the cross section and kinematics, which results 
in a net ~ 20% reduction in event rates at the current 
SK (and 10% reduction at a gadolinium-enhanced SK). 
While the SK limit could be weakened by these correc- 
tions [154| |. we adopt the published limit [34], since the 
net effect of the corrections is smaller than the statis- 
tical uncertainties on the number of background events 
at SK. Furthermore, corrections should also be applied 
to the source, which would reduce the neutrino opacity 
and lead to more high-energy neutrino escaping, partially 
canceling the effect of the corrections. One should also 
consider that there are uncertainties in approximating 
the neutrino spectrum, which will not be purely thermal. 

We show four DSNB spectra in the figure, as labeled. 
We also show the region already excluded by SK [34], EH 
by hatched shading. The exclusion region is drawn by 
scaling the 6 MeV spectra up so that the predicted event 
rate in the energy range 18-26 MeV is 2 events (22.5 
kton yr)" 1 . Since SK currently probes a narrow range of 
the exponential tail of the thermal neutrino spectra, the 
exclusion region has an almost temperature-independent 
shape. 

The SK limit already partially excludes the dominance 
of D e effective temperatures at or above 8 MeV. The 
6 MeV DSNB, which is typical of scenarios with neu- 
trino mixing, is within at most a factor 2 of the current 
SK limit. Furthermore, the 4 MeV spectrum, which re- 
flects our lower estimate for the effective temperature, 
lies within a factor 4 of the current SK limit. Effective 
temperatures in the relevant energy range lower than 4 
MeV are unlikely given the SN 1987A data and general 
considerations above. In Table HVl we show the predicted 
event numbers in the energy bin 18-26 MeV, which can 
be compared to the SK event limit, and which can be 
improved soon. 



TABLE IV: Integrated DSNB event rates in the positron en- 
ergy range 18-26 MeV for the current SK and 10-26 MeV 
for a gadolinium-enhanced SK. The uncertainties reflect the 
upper and lower i?ccsN, defined from the CSFH. 



v e spectrum 


events [(22.5 kton yr) -1 ] 




10 < E e /McV < 26 


18 < E e /MeV < 26 


8 MeV 


4.2 ± 1.4 


2.0 ±0.7 


6 MeV 


3.5 ±1.1 


1.3 ±0.4 


4 MeV 


1.8 ±0.5 


0.4 ±0.1 


SN 1987A 


1.7 ±0.5 


0.5 ±0.1 



C. DSNB future prospects 

In the current SK, DSNB D e are detected in a sin- 
gles positron search, D e p — > e + n, for which there are 
very large background rates [34[ • The largest background 
above ~ 18 MeV comes from the decay of invisible muons, 
i.e., non-relativistic muons produced from atmospheric 
and Dfj, [26| . The electrons and positrons produced in 
muon decays cannot be distinguished from those of the 
signal inverse beta decay reaction. At energies below ~ 
18 MeV, the decays of spallation products of cosmic-ray 
muons also contribute to the background (see Fig. [4j. 

With 0.2% dissolved gadolinium, the neutrons pro- 
duced in inverse beta decay can be identified with high 
efficiency, via the photons emitted upon neutron capture 
on gadolinium. This delayed signal allows tight tempo- 
ral and spatial coincidence for signal events, reducing the 
invisible muon background by a factor ~ 5 and remov- 
ing the spallation backgrounds in the range 10-18 MeV, 
opening it up to DSNB searches [45|,[46|]. This is shown in 
Fig. [5] where the relevant shadings shown in Fig. Q] have 
been removed. Below 10 MeV, reactor D e still overwhelm 
the DSNB signal 

Importantly, the enhancement allows a mostly rate- 
limited, rather than background-limited, DSNB search, 
so that the sensitivity improves linearly with exposure. 
Another important point is that a gadolinium enhance- 
ment could be applied to the existing SK detector [451 ] . 

The advantages of an enhanced SK are dramatic. At 
the lowered energy threshold, the event spectrum and 
uncertainties converge, so that the total range in predic- 
tions is only a factor ~ 2 (i.e., from the lowest predicted 
4 MeV curve to the highest predicted 6 MeV curve). The 
predicted total event rate in the energy range 10-26 MeV 
is almost triple that in SK with a high energy threshold, 
since the peak of the DSNB spectra can be probed; a 
gadolinium-enhanced SK is almost guaranteed to detect 
the DSNB. A non-detection would require novel stellar or 
neutrino physics, for example a change in the equation of 
state, extremely fast r otation, in visible neutrino decays 
on cosmological scales 156], or the effect of hypo- 

thetical particles on the emission model of supernovae. 
Any such explanations would have to show why the SN 
1987A detection results were atypical. 
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V. DISCUSSIONS AND CONCLUSIONS 

Neutrinos are the only probes (possibly apart from 
gravitational waves) of the central regions of core-collapse 
events, and their study and detection are strongly moti- 
vated by many areas of physics. In this paper we assess 
the uncertainties for DSNB searches and implications for 
stellar and neutrino physics. The results can largely be 
divided into three categories. 

A. On astrophysical inputs 

We start with an up-to-date compilation of the CSFH 
and cross-check it with other observables. While our aim 
is to constrain the astrophysical inputs for the DSNB 
(i.e., the CCSN rate), these new results are of value in 
their own right. 

• Consistency with CCSN rates: Using a collection of 
recently observed CCSN rates, we show that they 
are in agreement with predicted CCSN rates from 
the CSFH, up to z ~ 1. Importantly, this check is 
almost independent of the IMF. 

• Consistency with EBL: Using new EBL measure- 
ments, including TeV gamma-ray constraints, we 
show that predictions from the CSFH are in 
agreement with observations. The observed to- 
tal EBL is 73±26 nWm^sr 1 , while we predict 
78t^nWm- 2 sr" 1 for the BG IMF. 

• Consistency with stellar mass density: The choice 
of the BG IMF is independently supported by stud- 
ies of the stellar mass density, which show agree- 
ment in z < 0.7. 

In conclusion, our fiducial CSFH together with the 
BG IMF gives excellent agreement among all observa- 
tions considered, preparing the astrophysical inputs for 
the DSNB. In particular, we stress how the astrophysical 
inputs cannot be too low, in order to maintain the self- 
consistent picture of the birth, life, and death of stars, as 
illustrated by the cross-checks. 

B. On emission inputs 

Next we assess the range of neutrino emission expected 
from supernovae. Coupled to our cross-checked astro- 
physical inputs, this allows us to discuss how constraining 
the present limit on the DSNB is. Our best-determined 
results are shown in Fig. [5l 

• Neutrino emission: Ultimately, this is a quantity 
to be studied using supernova neutrino detections. 
We thus discuss the generic emission of 9 e from su- 
pernovae using simple arguments based on the en- 
ergetics of core collapse, neutrino mixing, and the 



observed SN 1987A neutrino burst. We also show 
the neutrino temperature depends weakly on the 
remnant parameters. We conclude that the rele- 
vant high-energy 9 e emission cannot be made too 
small. 

• Present constraints: The SK limit is already prob- 
ing interesting parameter regions. The dominance 
of an effective D e temperature above 8 MeV is con- 
strained. Similarly, a high rate of dark core col- 
lapses is also prohibited by an amount that depends 
on the assumed temperature. 

• Implications for stellar and neutrino physics: If an 
effective temperature is ruled out, this rules out 
neutrino mixing and initial neutrino temperatures 
that would lead to that effective temperature in 
the DSNB energy range. The correspondance be- 
tween the initial and effective neutrino emissions is 
dependent on the neutrino mixing scenario, which 
has been studied systematically by various authors. 

To conclude, from general and SN 1987A considera- 
tions, the neutrino emission cannot be too low. The cur- 
rent SK limit is already constraining interesting neutrino 
effective temperatures. Noting the SK limit was placed 
in 2003, more data and improved cuts will allow better 
sensitivity. 

C. On detection in SK 

With the DSNB inputs and their uncertainties checked, 
we discuss implications for DSNB detection in the future. 

• SK prospects: The DSNB is near detection in SK. 
The 6 MeV spectrum, typical of scenarios with neu- 
trino mixing, is within a factor ~ 2 of the current 
SK limit. The factor increases to ~ 4 for the lower 
4 MeV and SN 1987A reconstructed spectra. 

• SK with gadolinium: Intriguingly, the fluxes and 
uncertainties converge at the improved detection 
threshold (~ 10 MeV), so that predictions span an 
uncertainty of a factor ~ 2. The predicted event 
rate between 10-26 MeV is 1.2-5.6 events yr _1 , 
leaving no room for the DSNB to escape detection. 

• Future physics with SK: The effective temperature 
T Be contains physics concerning neutrino emission 
from the collapsed core and on neutrino mixing. 
Stellar and neutrino physics can be extracted by 
future comparisons between the measured and the- 
oretical Ty c . Future SK analysises should report 
results directly in terms of the time-integrated lu- 
minosity and effective temperature [46|, using the 
astrophysically-measured supernova rate. 

To conclude, while the current SK will continue prob- 
ing interesting physics, a gadolinium-enhanced SK is 
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almost guaranteed to detect DSNB events. A non- 
detection would require novel stellar or neutrino physics. 
Together with the decreasing astrophysical uncertainties, 
these results strongly support the imminent detection of 
the DSNB and solidify its important role in understand- 
ing supernova and neutrino physics. 
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